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Abstract 

Blending schemes based on circles provide smooth 'fair' interpolations 
between series of points. Here we demonstrate a simple, robust set of 
algorithms for performing circle blends for a range of cases. An arbitrary 
level of G-continuity can be achieved by simple alterations to the underly- 
ing parameterisation. Our method exploits the computational framework 
provided by conformal geometric algebra. This employs a five- dimensional 
representation of points in space, in contrast to the four-dimensional rep- 
resentation typically used in projective geometry. The advantage of the 
conformal scheme is that straight lines and circles are treated in a single, 
unified framework. As a further illustration of the power of the conformal 
framework, the basic idea is extended to the case of sphere blending to 
interpolate over a surface. 
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1 Introduction 

In a range of applications we often seek curves and surfaces that have an aesthet- 
ically pleasing 'roundedness' to them. One way to make this concept concrete 
is through looking for globally-optimised 'minimum variation curves' pp. The 
philosophy behind this idea is straightforward. We usually prefer curves that 
are close to circular over curves with sharp turns. This is particularly true when 
designing camera trajectories, where sudden changes in curvature can have a 
very disorienting effect. Circular paths are characterised by having constant 
curvature, so a natural idea in forming interpolations between control points 
is to find a curve that minimises the total change in curvature. The problem 
with such a strategy is that these curves can be extremely hard to compute. If 
one adopted a variational strategy, with endpoint conditions, the equations for 
the curve can be as high as fifth order and are even more difficult to treat than 
those of elasticity. Such equations can only be solved numerically and do not 
have straightforward, controllable, analytic solutions. The problem is even more 
acute if multiple control points are involved, as even numerical computation can 
be extremely difficult. 

A more straightforward, local scheme that provides smooth interpolations 
was introduced by Wenz [2] and later extended by Szilvasi-Nagy & Vendel 
The idea explored by these authors is to generate curves that are as close as 
possible to circles. Given four points Xq, . . . , ATa we construct the circles Ci 
through Xq, Xi and X2, and C2 through Xi, X2 and X3. The curve between 
Xi and X2 is then formed by smoothly interpolating between the two circles. 
This idea was further extended by Sequin & Yen and Sequin & Lee [HI, 
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who introduced an angle-based circle blending scheme. The angle-based scheme 
gives better results than the earlier, midpoint scheme, and we argue here that 
it is geometrically the 'correct' one. Sequin & Lee also showed how to achieve 
G'2-continuity (and higher order continuity, if desired), and demonstrated the 
value of angle-based blending for interpolating over the surface of a sphere. 

In this paper we further explore the geometry associated with circle blend- 
ing, following the methods developed by Sequin and his coworkers. The essential 
idea is that the natural way to transform between two circles is via a conformal 
transformation. Conformal transformations leave angles invariant, but can alter 
distances. Euclidean transformations are the subset of conformal transforma- 
tions that also leave distances invariant. Conformal transformations in a plane 
can take any three chosen points to any three image points. As such, they can 
transform a line or circle into any other circle. In this geometry, straight lines 
are examples of circles that pass through the point at infinity. By exploiting the 
features of conformal geometry, we can write robust code that treats (straight) 
lines and circles in a single, unified manner. This eliminates the need to check 
for special cases. Similarly, in three dimensions, planes and spheres are treated 
as examples of the same object. So a single routine can interpolate between 
points on a sphere, and will reduce to the planar case when four points happen 
to lie in a plane. 

To fully exploit the advantages of conformal geometry we work in the math- 
ematical framework of geometric algebra [HIE]. This algebra treats points, lines, 
circles, planes and spheres, and the transformations acting on them, in a uni- 
fied algebraic framework. A number of authors have argued for the advantages 
of the conformal geometric algebra framework for computer graphics applica- 
tions [71 IHl [Till The present work should be viewed in this context. We 
show how complex problems such as finding the conformal transformation be- 
tween a line and a circle reduce to simple, robust expressions in the geometric 
algebra framework. As a further application we show how the same framework 
naturally extends to sphere-blending over a surface. This suggests a new method 
of characterising surfaces that does not require the concept of swept curves. 

This paper starts with an introduction to conformal geometric algebra. This 
introduction is self-contained, but to keep its length down a number of concepts 
are introduced with a minimum of explanation. We then turn to the question 
of how to mathematically encode transformations between circles. We find the 
conformal transformation that achieves this and explore its properties. Some 
subtleties involving the orientation of the transformation are explained, and 
we demonstrate how they are easily resolved in the conformal framework. We 
then provide a series of examples of blended curves, and illustrate the effects 
of demanding higher-order G-continuity. We finish by introducing a method of 
sphere blending and discuss the potential of this idea for encoding surfaces. 

2 Conformal geometry and Euclidean space 

The starting point for our description of geometry is the conformal group. This 
marks a radical departure from conventional descriptions of Euclidean space 
based on projective geometry and homogeneous coordinates. The main advan- 
tage of basing the description in a conformal setting is that distance is encoded 
simply, making the geometry well suited to describing the real three-dimensional 
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world. 

Suppose we start with an n-dimensional Euclidean vector space K". The 
conformal group consists of the set of all transformations of R" that leave angles 
invariant. These include translations and rotations, so the conformal group 
includes the set of Euclidean transformations as a subgroup. The conformal 
group on M" has a natural representation in terms of rotations in a space two 
dimensions higher, with signature {n+l, 1). So, in the same way that projective 
transformations are linearised by working in a space one dimension higher than 
the Euclidean base space, conformal transformations are linearised in a space 
two dimensions higher. The conformal representation of points in Euclidean 
three-space consists of vectors in a 5-dimensional space. While this may appear 
to be an unnecessary abstraction, working in this fivc-dimcnsional space does 
bring a number of advantages. 

To exploit the conformal representation we need a standard representation 
for a Euclidean point in the five-dimensional conformal space. Given that we 
are occupying a space two dimensions higher, two constraints are required to 
specify a unique point. The first of these is that our underlying representation is 
homogeneous, so X and AX represent the same point in Euclidean space. The 
second constraint is that the vector X is null, 

X^^O. (1) 

(The existence of null vectors is guaranteed by the fact that the conformal 
vector space has mixed signature.) This is essentially the only further constraint 
that can be enforced which is consistent with homogeneity and invariant under 
orthogonal transformations in conformal space. Now suppose that ei, 62 and 
63 represent three vectors in the three-dimensional base space, and we add to 
these the vectors cq and 64. These satisfy 

el = -1, el = +1, (2) 

and all 5 vectors {gq, ... ,64} are orthogonal. From the two extra vectors we 
define the two null vectors n and n by 

n = 64 -I- Co, n = 64 — eo, n-n = 2. (3) 

(It is a straightforward exercise to confirm that these two vectors both have zero 
magnitude.) From these we need to chose a vector to represent the origin. This 
is conventionally taken as —■^n. The vector X can therefore be written as 

X ^ 2x -n + an (4) 

where x is the equivalent three-vector representation of the point in Euclidean 
space. The variable a must be chosen so that X is null, which fixes a — x^. We 
therefore arrive at the following representation of a point in conformal space: 

X = 2x + x^n - n. (5) 

This representation can be arrived at using more geometrical reasoning by con- 
sidering a stereographic projection [HJ |7]. It is immediately clear from equa- 
tion ((Sj that n represents the point at infinity. 
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Figure 1: The null cone. In a three dimensional space of signature (2, 1) the set 
of null vectors form a cone. To ensure a unique representation of points that 
maintains orientation, we restrict to the subsection with positive cq component. 
The set of vectors satisfying X ■ n — —2 defines the 'standard' representation of 
points, in this case defining the conformal representation of a one-dimensional 
space. 

The Euclidean coordinates of the point x are recovered from X via the 
homogeneous relation 

Xi = --V~^^ 1 = 1,2,3. (6) 
X -n 

This relation confirms that X and AX define the same Euclidean point. In 
principle, one could allow for A to be negative, but in practice this should be 
avoided. Restricting to positive A is is equivalent to stating that X must always 
have a positive cq component. The reason for this restriction can be seen in 
figure n The equation X^ = generates a cone structure, and we only want to 
employ one half of this to represent points. Maintaining a positive eo component 
throughout all algorithms enables us to keep track of orientations consistently. 

Given a point a, represented by the unnormalised vector A, we can place A 
in the standard form of equation (0) with the map 

One should be wary of employing this map when writing code, as the right- 
hand side is singular if A happens to be the point at infinity. As is the case for 
projective geometry, it is better practice to let the normalisation run free, and 
only use equation JSJ in the final stage to recover the coordinates. 

The power of the conformal representation starts to become clear when we 
consider the inner product of points. Suppose that X and Y are the conformal 
representations of the points x and y respectively, both in the standard form of 
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equation jSJ. Their inner product is 

XY = {x^n + 2x-n)- {y^n + 2y - n) 
= -2x^ - 2y^ + Ax-y 

^~2[x^yf. (8) 

This is the essential resuh that underpins the conformal approach to Euchdean 
geometry. The inner product in conformal space encodes the distance between 
points in Euclidean space. This is the reason why points are represented with 
null vectors — the distance between a point and itself is zero. 

Generalising equation ||SJ) to unnormalised vectors, the distance between 
points can be written 

\x-y\^^-2 ^-l . (9) 

Any transformation of X and Y that leaves this product invariant must represent 
a symmetry of Euclidean space. Because such a transformation must leave the 
inner product invariant, it must be an orthogonal transformation in conformal 
space. Such a transformation ensures that null vectors remain null vectors, and 
so continue to represent points. The transformation must also leave n invariant 
if the distance is to remain unchanged. We can now see that the Euclidean 
group is the subgroup of conformal transformations that leaves invariant the 
point at infinity. In what follows we will find applications for both Euclidean 
and more general conformal transformations. 



3 Conformal Geometric Algebra 

Geometric algebra was introduced by the nineteenth century mathematician 
W.K. Clifford, and has found many applications in physics [01 El- Often in 
this work the name Clifford algebra is used, but for applications in geometry it 
is becoming increasingly common to see Clifford's original name of geometric 
algebra to describe the field. A geometric algebra is constructed over a vector 
space with a given inner product. The geometric product of two vectors a and 
b is defined to be associative and distributive over addition, with the additional 
rule that the square of any vector is a scalar, 

aa = a^ e R. (10) 

If we write 

ab + ba= (a + b)^ - a'^ -b'^ (11) 

we see that the symmetric part of the geometric product of any two vectors is 
also a scalar. This defines the inner product, and we write 

a-b^ ^{ab + ba). (12) 

The remaining, antisymmetric part of the geometric product returns the outer 
or exterior product familiar from projective geometry. We write this as 

aAb= ^(ab-ba). (13) 

The totally antisymmetrised sum of geometric products of vectors defines the 
exterior product in the algebra. 
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The geometric product of two vectors can now be written 

ab = a-b + aAb. (14) 

Under the geometric product, orthogonal vectors anticommute and paraUel vec- 
tors commute. The product therefore encodes the basic geometric relationships 
between vectors. Now that we know how to multiply vectors together, it is 
straightforward to construct the entire geometric algebra of a given vector space. 
This is facilitated by introducing an orthonormal frame of vectors {e^}. These 
satisfy 

6i " Sj = Tjij , (15) 

where rjij is the metric tensor. For a space of signature (p, q), rjij is a diagonal 
matrix consisting of p +ls and q —Is. The space of interest to us here is the 
conformal space of signature (4, 1). A basis for this is provided by the vectors 
eo,...,e4, with eo having negative square. The algebra generated by these 
vectors has 32 terms in total, and is spanned by 

1 {ci} {ciAej} {ejACjAefc} {/ej / 
grade 1 2 3 4 5 

dimension 1 5 10 10 5 1. 

We refer to this algebra as Q{4,1). The term 'grade' is used to refer to the 
number of vectors in any exterior product. The dimensions of each graded 
subspace are given by the binomial coefficients. 

The highest grade term in ^7(4, 1) is called the pseudoscalar and is given the 
symbol /. This is defined by 

7 = 6061626364. (16) 

The pseudoscalar commutes with all elements in the algebra, a feature of odd- 
dimensional algebras, and the (4, 1) signature of the space implies that the 
pseudoscalar satisfies 

/2 = -1. (17) 

So, algebraically, / has the properties of a unit imaginary. But it also plays a 
definite geometric role, as multiplication by the pseudoscalar performs a duality 
transformation in conformal space. A matrix representation of ^7(4, 1) can be 
constructed in terms of 4 x 4 complex matrices. These can be found in the physics 
literature in the guise of the Dirac matrices [T3]. For practical applications 
this matrix representation has little value and one is better off coding up the 
algebraic rules explicitly. 

A general element of 5(4, 1) is called a multivector and can consist of a 
sum of terms all grades in the algebra. Arbitrary elements of C/(4, 1) can be 
added and multiplied together. A multivector that consists only of terms of 
a single grade is said to be homogeneous. The geometric product of a pair of 
homogeneous multivectors decomposes as follows: 

ArBs = {ArBs)\r—s\ + {^rBs)\r-s\ + 2 + ' ' ' + {ArBs)r+s- (18) 

The angle brackets {M)r are used to denote the projection onto the grade-r 
terms in M. The dot and wedge symbols are used to generalise the inner and 
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outer products to the lowest and highest grade terms in a general product: 

Aj--Bg = (y4ri3s)|i — s| 
ArhB,= {ArB,)r+s, (19) 

where and Bs are homogeneous multivectors of grade r and s respectively. 
One only ever uses the dot and wedge symbols when applied to homogeneous 
multivectors. Of particular importance are the inner and outer products with a 
vector. These satisfy 

a-Ar = \{aAr - {-ly ArO) 

a^Ar^\{aAr + {-lYAra). (20) 

Some straightforward algebra confirms that the outer product is associative . 

An important algebraic operation applied to multivectors is reversion. This 
plays an analogous role to transposition in matrix algebra. The reverse of a 
multivector M is denoted by M and is defined by reversing the order of all 
geometric products of vectors in M . An arbitrary multivector in 5(4, 1) can be 
written as 

M = a + a + B + T + Ib + Il3 (21) 

where a and f3 are scalars, a and b are vectors, i? is a bivector and T is a 
tri vector. The reverse of M is then given by 

M = a + a-B-T + Ib + I(3. (22) 

A conformal transformation in the Euclidean space MP is represented by an 
orthogonal transformation in 5(4,1). Of particular relevance here are special 
conformal transformations. These have determinant +1, and correspond to 
orientation-preserving transformations on R'^. In geometric algebra a special 
conformal transformation, applied to an arbitrary multivector A, can be written 
in the compact form 

A RAR. (23) 
Here R is an even-grade element (grades 0, 2 and 4) satisfying 

RR = 1. (24) 

Because our representation is homogeneous, this normalisation constraint can 
usually be relaxed to allow for RR to be an arbitrary positive scalar. The 
element R is called a rotor. As an example, the (unnormalised) rotor for the 
Euclidean translation between x and y can be written 

T^y ^ {n-Y)nX + {n-X)Yn, (25) 

where X and Y are the conformal equivalents of x and y respectively. Every 
rotor can be written in term of the exponential of a bivector as 

i? = ±cxp(-S/2). (26) 

The bivector B is the generator of the rotation, and is an element of the Lie 
algebra of the conformal group. 
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4 Geometric primitives in conformal algebra 



We have now established how points in R"^ are represented by null vectors in the 
conformal algebra ^/(4, 1). We now need to establish how higher dimensional 
objects are represented. Given a null vector Y, one view of the point this 
represents is as the solution space of the equation 

YAX = 0, X^ = 0. (27) 

where here X is treated as a vector variable. This is a similar idea to projective 
geometry, and generalises immediately to higher-grade objects. Given a mul- 
tivector Ar, the outer product of r distinct vectors, the geometric object that 
this represents is the solution space of 

ArAX = 0, X'^ = 0. (28) 

This pair of equations is clearly homogeneous in X . If a conformal transforma- 
tion is applied to Ar, taking it to RA^R, then the solution space transforms in 
exactly the same way (that is, X i-^- RXR). In this manner conformal transfor- 
mations are easily applied to higher-grade objects. 

After vectors, the next objects to consider are bivectors. Given a bivector 
B (formed from the outer product of two vectors), the solution space of equa- 
tion 1)28(1 depends on the sign of B^. If is negative there are no solutions, 
if = there is one solution, and if is positive there are two solutions. 
This final case corresponds to the situation where B = X f\Y , where X and Y 
are a pair of null vectors. Given a bivector of this form it is straightforward to 
recover X and K j^j, so a bivector can encode a pair of points. 

Next consider a trivector L. The main example of relevance is the trivector 
formed from the outer product of three points Xi, X2 and X3, so 

L ^ X1AX2AX3. (29) 

Any vector X that is a solution of LAX = must be a linear combination of 
Xi, X2 and X^. The fact that X is null and has arbitrary scale implies that 
in M'^ the solution space is one-dimensional. To see what this space is, consider 
the null vector C, where 

C = LnL. (30) 
If X satisfies X AL — 0, then X commutes with L. It follows that 

X-C 



X-nC-n C-n 



constant. (31) 



The solution set X therefore consists of points at an equal distance from some 
point C. It follows that the trivector L represents the circle through Xi, X2 and 
X3. If this circle also passes through the point at infinity, the line is straight. 
That is, the test that a line L is straight is that 

LAn = 0. (32) 

Similarly, the straight line through Xi and X2 is given by 

L = XiAX2An. (33) 



8 



Conformal algebra therefore treats lines and circles in a unified manner as trivec- 
tors in C/(4, 1). This is to be expected, as a conformal transformation can always 
transform a line into a circle. 

The radius p of the circle L is found from 

So again we see how a metric object (the radius) is recovered from a homoge- 
neous representation of a geometric entity (the circle L). The angle between 
two circles (or lines) is found from 

'""'^^wm- ^^^^ 

As expected, this is invariant under the full conformal group. 

Similar considerations apply to the 4-vector 5* defined by the four points 
Xi, . . . , Xi , 

S* = (36) 

The object defined by S is the unique sphere through the four points (which 
cannot all lie on a line or circle if 5* 7^ 0). The centre of the sphere is given by 
SnS, and if the centre lies at infinity the sphere reduces to a plane. As with the 
case of a straight line, the plane defined by the three points Xi, X2 and X^ is 

p = Xl^X2^Xz^n. (37) 

Lines, circles, planes and spheres can all be intersected in a straightforward 
manner in conformal geometric algebra. The two cases of greatest significance 
are those of a circle and a sphere, and of two spheres. Treating the latter first, 
any two spheres intersect in a circle (provided they touch). This reduces to a 
line if both spheres are in fact planar. The intersection L is found simply by 

L={ISi)-S2, (38) 

which is a trivector, as required. Similarly, the intersection of a circle C and a 
sphere S is given by 

B^{IC)-S. (39) 

This is a bivector, as a circle can intersect a sphere or plane in zero, one, or 
two points. In the case where C is a straight line and S" is a flat plane, one 
of the points of intersection is at infinity. As we proceed we will see how the 
orientations implicit in lines and surfaces are neatly encoded in these intersection 
formulae. 



5 Circle blending 

Suppose we are given a series of rt + 1 points xq, . . . , a;„ and we seek a smooth 
curve through all of these points. The idea behind the circle blending scheme is 
to find a curve that is as close as possible to a circle and which passes through 
all of the points. A typical setup is illustrated in figure |21 Consider the four 
points Xi^i . . .Xij^2- We form the circle Ci through x^-i, Xi and Xi+i, and the 
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Figure 2: Circles splines in the plane. The curve through each intermediate 
point is generated by blending between the two circles at each end, producing 
a smooth curve with limited change in its curvature from point to point. The 
curve shown has G2 continuity, using the scheme described in section [7| 



circle C^+i through x^, x^+i and In the region between xi and Xi^i we 

need a curve that blends smoothly between the circles Cj and Ci+i. 

A range of schemes exist for interpolating between circles, but one is natu- 
rally picked out from the viewpoint of conformal geometry. We are presented 
with two circles Ci and C2, both of which share two points in common: 

Xl^X2^Xy,. (40) 

Given the (conformal) points Xi, X2^ and the circles Ci and C2, there are no fur- 
ther objects to specify. First we need to define family of circles between Ci and 
C2. This is easily achieved by finding the transformation between them. The 
transformation must map between two circles, so is a conformal transformation. 
The only bivector generator for such a transformation is the antisymmetrised 
product between Ci and C2. A transformation governed by this generator is 
simply a rotation in Q{A, 1). If we normalise the circles by defining 

then the interpolated circle is simply 

C12W = (sin((l - A)^)^! + ^xn{\e)C2) • (42) 

Here 9 is the angle between the circles, defined by 

cos(6i) = C'i-C'2. (43) 

(The case of 6* > tt will be addressed shortly.) This method of interpolating 
between circles recovers the angle-blending scheme of Sequin & Yen thus 
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Figure 3: Interpolating between a pair of circles. The two outer circles share two 
points in common, and between them specify a sphere. AU of the intermediate 
circles also lie on this sphere. 

providing a firm geometric reason for preferring that scheme. At this point 
there is no reason for Ci and C2 to lie in the same plane. If they lie on different 
planes then the four points X^, . . . , ATa define a sphere. Since the interpolated 
circles can be written 

Ci2(A) = -i^Xi (^sin((l - X)e) ^ + sin(A0)^) , (44) 

we see that all points on these circles are combinations of Xq . . . X^. It follows 
that all of the intermediate circles also lie of the sphere defined by Xq, . . . , X^, 
which is a desirable property for a range of applications 5 . The basic interpo- 
lation scheme is illustrated in figure 01 

Now that we have a straightforward means of encoding the circle blends, we 
need to parameterise the actual trajectory between the points. A simple means 
for achieving this is to start with the straight line between Xi and X2. This 
path is parameterised by 

Yi2{X) = -(1 - X)X2-nXi - \Xi-nX2 + X{l - X)Xi-X2n, (45) 

and the line itself is described by 

ii2 = ATiAXsAn. (46) 

The rotor that transforms between the line L12 and the circle Ci2{X) is simply 

i?i2(A) = l + Ci2(A)ii2. (47) 

So once the (normalised) blended circle 6*12 (A) is found, the path itself is given 
by transforming from the straight line to the circle. The path between Xi and 
Ar2 is therefore simply 

Xi2(A) - i?i2(A)ri2(A)i?i2(A). (48) 

This is extremely simple to code up. We never have to calculate the centre of a 
circle, and the algorithm deals equally easily with straight lines or circles. The 
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Figure 4: Circles splines on a sphere. The curve is generated using the same 
algorithm as for the planar case. The algorithm ensures that all intermediate 
points in the middle section lie on the sphere specified by the four control points. 
In the example, all five control points lie on the same sphere, so the curve sits 
smoothly on the sphere. 

algorithm keeps all intermediate points on the sphere defined by each group of 
four points. If all of the base points lie on the same sphere, then we compute 
smooth curves over this sphere. An example of this is shown in figure 0] 

6 Problems and special cases 

Provided the sequence of any four points are unique, the circles to be blended 
and the intermediate straight line are all well defined. The only aspects of the 
algorithm that can be problematic are concerned with the definition of the angle 
6 in equation (|42|l . There are two cases over which care must be taken, both 
of which can be treated fairly easily. The first obvious problem is that the 
definition of the blended segment breaks down if sin(0) — 0, which occurs when 
cos{9) — +1 and when cos(0) = —1. The former case corresponds to 6* = and 
implies that the two curves are the same. This is easily caught by setting a 
threshold value for 9 below which eqn 142() is replaced by the small-f? limit of 



The case oi 9 — ir is more complicated, and corresponds to a somewhat 
pathological example. This occurs when the two circles only differ in their 
orientation, so C2 = —Ci- A sample configuration is shown in figureEl AH four 
points lie in a plane, so the interpolated curve must also lie in this plane. It is 
possible to define an interpolation scheme for this case. If we let P denote the 
plane containing the circles, then the generator of the transformation must act 
in the P plane, leaving Xi and X2 invariant. The bivector generator for the 
transformation is therefore 



Ci2(A) = (1- A)Ci +AC2. 



(49) 



B12 = {XiAX2)P 



(50) 
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Figure 5: Four points on a circle. The four points lie on a circle, but the two 
circles we interpolate between are equal and opposite. An interpolation scheme 
in the plane is quite possible, but the case is pathological. If one point moves 
slightly off the plane, the interpolated curve will jump to a spherical case. 

This produces the type of interpolated curve shown in figure |S1 The problem 
here is that this case is not smoothly connected to the general case. Imagine 
moving one of the control points in figure El slightly out of the page. The four 
points then define a sphere, and the interpolated curve will lie on this sphere. 
This is a totally different curve to the planar case. If we are only interested in 
planar plots, then everything is well-behaved and the above case can be easily 
dealt with. But if we are interested in curves in three-dimensional space, this 
case is best avoided altogether with the addition of further control points. 

The case of oppositely oriented circles relates to the second key issue, which 
is what happens if the circles are greater than 180° apart? The arccos function 
will not return the correct angle, and the interpolation scheme will effectively 
blend the wrong way. This problem has an elegant solution within the conformal 
framework. The question reduces to finding the correct midpoint circle Ci^2- 
This is given by 



The case of 9 = n corresponds to equal and opposite circles, which is the one 
awkward case we have to avoid. All we need is a test to decide which is the 
correct mid-circle for the geometry, and the arccos function can be used to 
determine 9/2 unambiguously from 



To find a suitable test we consider the plane of points equidistant from Xi and 
X2- This plane is defined by 




9 <7r 
9>Tr 



(51) 



cos(6i/2) = Ci-Ci,2. 



(52) 



n = I{Xi/\X2)-n. 



(53) 
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Figure 6: Intersection of a plane and a circle. Both the plane and the circle are 
oriented objects, and the intersection bivector defined by equation H54(l returns 
XfAX,. 

All planes are oriented, so the sign here is important. This choice corresponds 
to the normal to the plane (defined by the right-hand rule) pointing in the 
direction from Xi to X2- Now suppose that we intersect a circle C with the 
plane H. The result is two points encoded in the bivector _B, where 

s = /(cn + nc). (54) 

The bivector B contains the points in the order 

B^XfAXi, (55) 

where ^/ is the point where the circle intersects the plane from the negative to 
positive side, and Xi is the opposite point (see figure IHl- 

Now suppose that Mi and M2 are the midpoints of the circle segments Ci 
and C2 respectively. Both lie in the plane 11. The correct mid-circle should have 
its Xf intersection point closer to the midpoint of Mi and M2 than its Xi point. 
If this is not the case, then the circle has the wrong orientation. We therefore 
only need form the quantity 

a= {(Mi-nM2 + M2-nMi)An)-B (56) 

where B is given by equation 154|l with C defined by 

C^Ci + C2. (57) 

The test is now 

Ci.= |^^ + ^^, (58) 
■ \-(Ci+C2), a>0 

This conditional statement covers all special cases (where various circles are 
degenerate or reduce to straight lines) and the border cases of a = is the 
single degenerate case mentioned above. With the mid-circle found, the angle 9 
is computed, and the interpolation scheme of equation H42|l can be followed. 

A further problem with the blending scheme is highlighted in figure |7| The 
geometry of the control points ensures that an unwanted inversion arises in the 
middle of the plot. This is because the gradient at each control point is fixed 
entirely by the two points on either side (see section [71), and so in this case is 
forced to be flat, as illustrated. The problem can be removed in a number of 
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Figure 7: An unwanted inversion. The blending scheme does not always min- 
imise curvature, and can lead to additional inversions in certain cases, as shown 
on the left-hand side. This is because the gradient at any point is determined 
entirely by the two points directly next to it. The problem can be resolved by in- 
troducing additional midpoints in the problem area and re-blending, producing 
the smooth curve shown on the right. 

ways. One of the simplest is to define a series of midpoints, and then repeat the 
blending algorithm. The result of a single iteration of this type is also shown 
in figure Of course, one can go further and continue to recursively define 
midpoints, which will also generate a smooth blend. This latter approach has 
the attractive feature that it removes the need for any trigonometric evaluations, 
as the mid-circle is found easily using the method described above. 

7 Continuity 

An important issue in any interpolation scheme is the order of continuity at the 
control points. In this section we will only consider G-continuity. To achieve 
C-continuity some reparameterisation will have to occur on each spline. For the 
case of a camera fly-by, this can be achieved by specifying the desired velocity 
along the curve. 

Our blending scheme is defined by equation iPHji . and before we can differ- 
entiate this we need a result for the derivative of a rotor. The definition of a 
(normalised) rotor is sufficient to prove that we can always write 

^ - -\RB (59) 

where B is a bivector, which in general will also be a function of A. (The factor 
of —1/2 is a useful convention). To find the tangent vector to the (conformal) 
curve X{X) we form the line 

T = X^X'^n (60) 

where the dash denotes the derivative with respect to A. To see how this behaves 
we return to equation H48|) and form the derivative at A = 0. For convenience 
we will assume that the rotor in equation H47|l has been normalised, though this 
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is unimportant when forming the tangent vector. In evaluating the derivatives 
we can use the resuhs that 

R{\)XiR{X) = Xi, R{X)X2R{\) = X2. (61) 

Differentiating these, we see that 

Xi-B{X)=X2-B{X) = 0, (62) 

and these also hold for derivatives of B with respect to A. 

On differentiating H47|l we find that (dropping the subscripts on R and Y) 

X' = R{Y-B + X2-nXi ~Xi-nX2 + (1 - 2X)Xi-X2n)R. (63) 

It follows that 

X'\^^a = X2-nXi-Xi-nX2 + Xi-X2RonRo, (64) 

where Rq = R{X = 0). The tangent vector at Xi is therefore given by (up to a 
scale factor) 

Ti = Xi-nX2AXiAn- Xi-X2 (RonRa)AXiAn 
= -(i?o((XiAX2An)-Xi)i?o)An 

= -{Ci-Xi)An. (65) 

But this is precisely the tangent vector to the circle Ci at that point Xi. So 
the tangent at each of the control points is defined by the circle through the 
control point and its two adjoining points. The tangent vectors at a connection 
are therefore continuous, which ensures Gi continuity. 

Next we need to consider G2 continuity. For this we need the circle 

a = xax'ax" 

^ r{ya{Y' + Y-B)a{Y" + 2Y' B + Y-B')^R. (66) 
Evaluating this at A = we obtain 

= RQ{-X2-nXi)A{X2-nXi- Xi-nX2 + Xi-X2n)A 

{-2Xi-X2n + 2Xi-X2n-B{Q))Ra 
= -2{Xi-nX2-nXi-X2)C\ 

-2{X2-nXi-X2)RoXiA{-Xi-nX2 + Xi-X2n)A{n-B{Qj)Ro. (67) 

This contains a term in Ci, and an additional term controlled by the bivector 
B (evaluated at A = 0). The first term is the desired one, as it ensures that 
the radius of curvature at the control point is defined by the circle through it. 
The second term is not wanted, and implies that the simple scheme defined 
by equation H42I) does not guarantee G2 continuity. This was first pointed out 
by Sequin & Lee 0. To provide G2 continuity we need only ensure that the 
derivative of the interpolating rotor R vanishes at the endpoints. This is simply 
achieved by replacing equation (|^ with 

C'i2(A) - fsin((l - 3A2 + 2X'')e)Ci + sin((3A2 - 2X'')e)C2) • (68) 

sm 6^ v ^ ' / 
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Figure 8: Comparison between Gi and G2 continuity. The Gi blend is shown 
on the left, and the G2 blend on the right. The G2 blend hugs the control points 
more closely than the Gi case, and for most cases the G2 curves are the more 
aesthetically satisfying. 

This blending scheme does now have G2 continuity, and is the scheme use to 
produce the figures in section [S] In figure |S1 we show a comparison between the 
two schemes. As expected, the G2 scheme smoothes out some of the curvature 
around the control points. One can extend this idea to obtain any desired order 
of continuity. This is achieved by replacing the polynomial in equation H68|l by 
a higher-order blend. 

8 Sphere blending 

One can straightforwardly apply the ideas developed here to swept surfaces 
using a version of the scheme described by Szilvasi-Nagy & Vendel [H]. But in 
this section we aim to explore an alternative idea, based on sphere blending. A 
sphere is described as a 4-vector in the conformal geometric algebra, and these 
can be transformed and interpolated in a similar manner to circles. 

As a simple example, suppose that the points Xi, X2 and define the 
vertices of a triangle. To each corner we attach a sphere, S'l, 5*2 and S3. Each 
sphere passes through the three vertices of the triangle, and a fourth point which 
can be viewed as a control point. That is, we can write 



These assume that none of Ai, A2 and A3 lie on the circle defined by Xi, X2 
and X3. We next need to define a blend over the surface. For this we introduce 
the barycentric coordinates (A, /i, ^'), subject to 



Si 
S2 
S3 



AiAXiAX2AX3 
yl2AXiAX2AX3 
A3AX1AX2AX3. 



(69) 



0<A,^,z^<l, A + /.t + i^ = l. 



(70) 
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We let F(A,/i) denote the conformal representation of the point in the triangle 
corresponding to the barycentric coordinates (A, /i, 1 — A — /i). Similarly, we 
define a (linear) sphere blend by 

S{X,fi) = XSi+ 1182 + 1^83. (71) 

We could employ a trigonometric blending scheme over the (abstract) sphere 
defined by the unit 4- vectors Si , S2 and 5*1 , but this raises a number of compli- 
cations. There is no single straightforward generalisation of barycentric coordi- 
nates over a spherical triangle, and each alternative scheme has its own merits 
and drawbacks 114|. Here we have adopted the simplest, linear blending scheme, 
which generates interesting surfaces. 

Now that we have the sphere and the point on the triangle defined, all that 
remains is to define the conformal transformation from the point to the blended 
surface. First we define the plane P through the three base points, 

P = XiAXaAXgAn. (72) 

Next we define the rotor R for the conformal transformation between the plane 
and the blended sphere, 

R{X,ti) = l-S{X,^i)P. (73) 

Finally, the surface itself is defined by the points X{X, /i), where 

X{X, ii) = i?(A, Ai)r(A, ii)R{X, ii). (74) 

A typical surface defined in this manner is shown in figure |51 The result is 
aesthetically quite pleasing, yielding a smooth blend free of sudden changes in 
curvature. 

Much work remains to extend this idea to a complete framework for defining 
blends over surfaces. The control points need to be chosen so as to reflect the 
geometry around the triangle, and continuity between blended spheres may be 
hard to achieve. Here we hope to have demonstrated that such an approach 
may be feasible. 

9 Summary 

Smooth splines between control points can be defined in terms of blended cir- 
cles. The splines pass through all points, and have no extra control points. 
The natural geometric framework for handling circles is conformal geometry. 
This geometry is encapsulated in a simple, unified manner by employing the 
geometric algebra of a five-dimensional space. All transformations are easily 
defined, and algorithms can be written in such a way as to minimise problems 
with special cases. Similar ideas can be applied to sphere blends over a surface, 
and in future work we will explore the potential of this idea for design and for 
encoding surface data. 
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Figure 9: A blended surface. A sphere is attactied to eacli of tlie vertices of a 
triangle. The surface is defined by a linear blend of the spheres. 
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